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Abstract. In the standard flat cosmological constant (A) cold dark matter (CDM) cosmology, a model of two populations of 
lens halos for strong gravitational lensing can reproduce the results of the Jodrell-Bank VLA Astrometric Survey (JVAS) and 
the Cosmic Lens All-Sky Survey (CLASS) radio survey. In such a model, lensing probabilities are sensitive to three parameters: 
the concentration parameter c\, the cooling mass scale M c and the value of the CDM power spectrum normalization parameter 
<x 8 . The value ranges of these parameters are constrained by various observations. However, we found that predicted lensing 
probabilities are also quite sensitive to the flux density (brightness) ratio q r of the multiple lensing images, which has been a 
very important selection criterion of a sample in any lensing survey experiments. We re-examine the mentioned above model 
by considering the flux ratio and galactic central Super Massive Black Holes (SMBHs), in flat, low-density cosmological 
models with different cosmic equations of state co, and find that the predicted lensing probabilities without considering q r 
are over-estimated. A low value of q r can be compensated for by raising the cooling mass scale M c in fitting the predicted 
lensing probabilities to JVAS/CLASS observations. To determine the cosmic equation of state co, the uncertainty in M c must be 
resolved. The effects of SMBHs cannot be detected by strong gravitational lensing methods when q r < 10. 

Key words, cosmology: theory - gravitational lensing. 



1. Introduction 

Gravitational lensing provides us with a powerful probe of 
mass distribution of the universe. By comparing the lensing 
probabilities predicted by various cosmological models and the 
density profile of lenses with observations, we are able to test 
the mass distribution of dark matter halos and, in particular, the 
inner density slope because the Jodrell-Bank VLA Astrometric 
Survey (JVAS) and the Cosmic Lens All-Sky Survey (CLASS; 
Browne et al. |2U0Tj| Helbig lSUUUl Browne et al. 2002, Myers et 
al. 2002 ) has provided us with observed lensing probabilities at 
small image separations (0.3" < A9 < 3"). 

In a standard flat ACDM cosmology, it is believed that a 
model of two populations of lens halos for strong gravitational 
lensing can reproduce the results of the JVAS/CLASS radio 
survey (Porciani & Madau l2000l Kochanek & White I200T1 
Keeton l200Tl Keeton & Madau l200ll Sarbu, Rusin & Ma l200Tl 
Li & Ostriker l200l 



Oguri I20021 Oguri I20031 Oguri, Lee 
& Suto 2003 1. The two populations of lens halos are distin- 
guished by introducing an abrupt change in the structure of 
the objects at the cooling mass scale M c dividing galaxies and 
clusters. Some authors used Singular Isothermal Sphere (SIS) 
models below M c and Navarro-Frenk- White (NFW) models 
above M c and could fit the observed image separation distri- 
bution of JVAS/CLASS (Li & Ostriker 1277011 Sarbu, Rusin & 
Ma 2001 1. Based on the fact that the image separation distribu- 
tion of lenses below 1 depends sensitively on both the inner 
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mass profile of galactic halos and the faint end slope of the 
mass and luminosity functions, Ma ( 2003 ) compared the tradi- 
tional approach that models the lenses as SIS and the Schechter 
luminosity function with a dark matter based approach that 
models the lenses with a certain halo mass profile and the 
Press-Chechter mass function. Constraints on the inner total 
mass profiles of halos are investigated by requiring the two ap- 
proaches to give consistent predictions. Li & Ostriker |2003l 
further proposed a model of three populations of halos as lenses 
distinguished by the halos mass to calculate the lensing proba- 
bility of image separation and time delay. 

In this paper, we revisit the model of the two populations 
of lens halos to calculate the lensing probabilities in a flat 
quintessence cold dark matter (QCDM) cosmology with dif- 
ferent cosmic equations of state co. The focus here is on the 
flux density ratio q r and cooling mass scale M c , and their ef- 
fects on the estimate of the equation of state of dark energy 
co. In spite of the fact that low mass lensing halos (galactic 
size or less) cannot all be SIS (Ma l^UOH Li & Ostriker |2"U0H 
Benson et al. 2002 Somerville & Primack 1999; Kaufmann et 
al. Qj93 ), we still model the galactic lens halos as SIS in our 
calculations based on the following considerations. On the one 
hand, evidence based on stellar dynamics of elliptical galaxies 
(e.g., Rix et al. 119971 Romanowsky & Kochanek 1999; Treu 
& Koopmans 2002 1, modeling of lensed systems (e.g., Cohn et 
al. 2001 1 and flux ratios of multiple images (Rusin & Ma l200li 
Rusin 2002 1 all give an inner profile for lensing galaxies that 
is consistent with SIS. On the other hand, since the SIS lensing 
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cross section is several orders of magnitude higher than that of 
the NFW, the lensing probabilities depend strongly on the cool- 
ing mass scale M c because a larger M c (and a smaller M c at the 
lower mass end of SIS lens halos, see Li & Ostriker 2003J al- 
lows more halos to be modeled as SIS. So the importance of the 
subdivision for low mass lens halos is lessened unless the large 
uncertainty of M c is sufficiently reduced. We thus ignore such 
a subdivision in this paper, and focus on the roles played by the 
flux density ratio q r and cooling mass scale M c in constraining 
the equation of state to. 

In addition to q r and M c , lensing probabilities estimated by 
the NFW halo are also sensitive to the concentration param- 
eter c\. We use the mass-dependent fitting formula given by 
Bullock et al. (2001 1. For a given halo mass and redshift, there 
is a scatter in c\ that is consistent with a log-normal distribu- 
tion with standard deviation tr c = A(logc) 0.18 (Jin^2000 
Bullock et al. 2001 ). Taking into account the scatter in c\ by av- 
eraging the lensing probability with the lognormal distribution 
will increase the probabilities considerably at larger image sep- 
arations and only slightly at smaller separations (Chen 2003b I. 
Since the well-defined sample of JVAS/CLASS is limited to 
small image separations (0.3" < AO < 3"), we thus ignore the 
scatter in c\ and our conclusions will not be affected. Another 
major uncertainty in predicting lensing probability arises from 
the considerable uncertainty regarding the value of the CDM 
power spectrum normalization parameter cr 8 (Chen 2003b i; 
we adopt <x 8 = 0.95, suggested by the abundance of massive 
clusters at redshifts z = 0.5 ~ 0.8 (Bahcall & Bode [2003) 
which gives <x 8 = 0.98 + 0.1. Lensing probability increases 
quickly with the source redshift z s \ since its distribution in the 
JVAS/CLASS survey is still poorly understood, we use the es- 
timated mean value of < z s >= 1-27 (Marlow et al. 2000 1. 

On the other hand, it should be pointed out that, like most 
authors, we use the spherical lensing profiles. Non-spherical 
lensing profiles do not change significantly the separation of 
observed multiple images, unless the profile deviates signifi- 
cantly from the spherical one. However, as pointed out by many 
authors (e.g. in a cluster environment, this has been discussed 
by Meneghetti et al. 12003 a 2003b, for the importance of sub- 
structures see Mao & Schneider 1998), the elliptical lenses and 
their substructures do have significantly larger cross sections. 
A more realistic lensing model should include these effects, 
which will be discussed in another paper. 

So in this paper we use the model of two populations of 
lens halos mentioned above to investigate the effect of image 
flux density ratio q x on lensing probabilities and find that the 
predicted lensing probabilities without considering q r are over- 
estimated. In flat, low-density cosmological models with dif- 
ferent cosmic equations of state to, we show that when the flux 
density ratio q x changes from 10 to 10 4 , the corresponding fit 
value of the cosmic equation of state will change from to = -1 
tow = -1/2. Also, in our flux-limited statistics of strong lens- 
ing image separations, the contributions from galactic central 
SMBHs can be safely ignored. 

The paper is organized as follows: Sect.|2]provides a brief 
description of the cosmology model used in this paper, Sect. [3] 
presents the predicted integral lensing probabilities, Sect. |4] 
gives our results, the discussion and comparisons with previ- 



ous work. Appendix A is devoted to a detailed deduction of 
lensing equations for galactic halos. 

2. Cosmological model and mass function 

In this section, we describe the cosmological model and dark 
halo mass function. More and more evidence shows that the 
Universe at present is dominated by a smooth component with 
negative pressure, the so called dark energy. So far, two kinds 
of dark energy have been proposed, quintessence and the so- 
called Chaplygin gas (e.g., Bento, Bertolami & Sen [2002 
2003|. We consider only the quintessence in this paper. It is 
assumed that quintessence offers an alternative to the cosmo- 
logical constant as the missing energy in a spatially flat uni- 
verse with a subcritical matter density Q m (Caldwell, Dave & 
Steinhardt fl998t . In this paper, we study spatially flat QCDM 
models in which the cold dark matter and quintessence-field 
make up the critical density (i.e., Q m + Q.q = 1). Throughout 
the paper, we choose the most generally accepted values of 
the parameters for flat, QCDM cosmology, in which, with the 
usual symbols, the matter density parameter, dark energy den- 
sity parameter and Hubble constant are respectively: Q m = 0.3, 
Q.q = 0.7, h = 0.75. Four negative values of to in equation 
of state pq = u>pq, with to = -1 (cosmological constant), 
to = -2/3, to = -1/2 and to — -1/3, are chosen to see their ef- 
fects on lensing probabilities. We use the conventional form to 
express the redshift z-dependent linear power spectrum for the 
matter density perturbation in a QCDM cosmology established 
by Ma et al. GWM 



P(k,z) = A Q k n Tl(k,z) 



8q(z) 



*q(0)(1+z) 



(1) 



The normalization Aq can be determined by cr 8 for different 
values of to 



1 f°° 

— r P(k,0)W 2 (kr s )k 2 dk, 
2n z Jo 



(2) 



in which W(x) = 3[sin(jt)/x 3 - cos(x)/jt 2 ] is the Fourier trans- 
formation of a top-hat window function, and r 8 = 8/i _1 Mpc. 
We choose the spectral index of the primordial adiabatic den- 
sity perturbations to be n = 1 . Tq is the transfer function for 
QCDM models, which is related to T\, the transfer function 
for ACDM models, with the relative transfer function 



Tqa 



Tq a + aq 2 



1 + aq 



2 ' 



q = 



(3) 



where Fq is the shape parameter in QCDM models and a is 
a coefficient that quantifies the relative amplitude of the mat- 
ter density field on large and small length scales, and is well 
approximated by 



a = (-to) s , 

s = (0.012 - 0.036a. - 0.017/w)[l - Q. m (z)] 
+(0.098 + 0.029w - 0.085/w)lnO m (z), 



(4) 



where Q m (z) = Q m /[Q m + (1 - Q m )(l + z) 3< "] is the matter den- 
sity parameter. We use the fitting formula for ACDM transfer 
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function given by Eisenstein & Hu ( 1999 1: T A = L/(L + Cq 2 eff ), 
with L = ln(e + 1.84^ eff ), q eS = £/(Q m /! 2 Mpc~') and C = 
14.4 + 325/(1 +60.5^"). 

Similarly, the linear growth suppression factor of the den- 
sity field gQ for QCDM in equation ([0 is related to that for 
ACDM g A with 

8Q , 

t = -(0.255 + 0.305w + 0.0027 /<j)[1 - Q m (z)] 

-(0.366 + 0.266w - 0.07/w)lnQ m (z), (5) 

where the empirical fit of gA is given by gA = 
2.5Q m (z)/{Q m (z) 4/7 - 1 + Q m (z) + [1 + Q m (z)/2][1 + (1 - 
^m(z))/70]} (Carroll, Press & Turner fT992l . 

We know that most of the consequences of quintessence 
follow from its effect on evolution of the expansion rate of the 
Universe, which in turn affects the growth of density pertur- 
bations (as described above) and the cosmological distances. 
From the first Friedmann equation dipper 1 ) = -pQda 3 (a = 
1/(1 + z) is the scale factor of the Universe), the expansion rate 
(Hubble constant) for a flat Universe and then the proper dis- 
tance and the angular-diameter distance can be calculated in 
our QCDM models (Huterer & turner 2001 1, which are needed 
in predicting lensing probabilities. 

The physical number density h(M, z) of virialized dark ha- 
los of masses between M and M+dM is related to the comoving 
number density n(M,z) by h(M, z) = n(M, z)(l + z) 3 , the latter 
originally given by Press & Schechter ( 1974b . and the improved 
version is n(M,z)dM - (p(0)/M)f(M,z)dM, where p(0) is the 
current mean mass density of the universe, and f{M, z) is the 
mass function for which we use the expression given by Jenkins 
et al. ll200TV 

3. Lensing probabilities 

In this section, we first give a brief description of lensing equa- 
tions for galaxies and clusters of galaxies, then present the pre- 
dicted integral lensing probabilities. According to the model 
of two populations of halos, cluster-size halos are modeled as 
NFW profile: pnfw = Ps>~l /[r(r+r s ) 2 ] , where p s and r s are con- 
stants. We can define the mass of a halo to be the mass within 
the virial radius of the halo r v ; r : Mum = 47rp s r 3 /(ci), where 
/(c0 = ln(l + ci) - ci/(l + ci), and c x = r vir /r s = 9(1 + 
zy l (M/l.5 x lO^Zr'Af©) -0 ' 13 is the concentration parameter, 
for which we have used the fitting formula given by Bullock 
et al. (2001 1. The lensing equation for NFW lenses is as usual 
y = x - n s g(x)lx (Li & Ostriker 120021 Chen I2003al I2003bi . 
where y - [y|, i\ = yD^/D^ is the position vector in the 
source plane, in which and are angular-diameter dis- 
tances from the observer to the source and to the lens respec- 
tively, x = and x = £/r s , £ is the position vector in the lens 
plane. Since the surface mass density is circularly symmetric, 
we can extend both x and y to their opposite values in our actual 
calculations for convenience. The parameter p s = 4p s r s /£ cr is 
x independent, in which S cr = (c 2 /4ttG)(D£/D^D£ s ) is the 
critical surface mass density, with c the speed of light, G the 
gravitational constant and D^ s the angular-diameter distance 



from the lens to the source. The function g(x) has an analytical 
expression originally given by Bartelmann ( 1996). 

Since the observational evidence presented so far suggests 
the ubiquity of black holes in the nuclei of all bright galaxies 
regardless of their activity (Magorrian et al. Q998 , Ferrarese 
& Merritt l2000l Ravindranath et al. 1200 II Merritt & Ferrarese 
200 fal l2001bl 1200 lcl Wandel l2002l Sarzi et al. 120021 we add 
a central supermassive black hole (SMBH) as a point mass to 
each SIS modeled galactic lens halo in its center(Keeton 2002 1 
Chen l2003all2003b> . It is well known that the finite flux den- 
sity ratio will definitely reduce the value of the lensing cross 
section(Schneider, Ehlers & Falco 1992). In the model of one 
population of halos (NFW) combined with each galactic halo 
a central point mass, the lensing probabilities are shown to be 
sensitive to the flux density ratio (Chen 2003a 2003b). One 
may argue that this is the case only because we treat the galac- 
tic bulge as a point mass (some fraction of the bulge mass), 
since a point mass lens will always produce a larger flux density 
ratio q x of the multiple images. So it would be interesting to in- 
vestigate the model of two populations of halos to see whether 
or not the lensing probabilities will be sensitive to q T , We will 
show later that the answer is yes, and the contribution from an 
individual central galactic SMBH is detectable only if q T > 10 3 . 
When q T < 10 (for the well-defined sample of JVAS/CLASS), 
the effect of SMBH can be safely ignored. 

So, for galaxy-size lenses, we use an SIS+SMBH model. 
The surface mass density is 

= M.6 2 (£) + (6) 

where M. is the mass of a SMBH, £ is the position vector in 
the lens plane, £ 2 (£) is the two dimensional Dirac-delta func- 
tion, and o" v is the velocity dispersion. Choosing certain length 
scales, the lensing equation can be written as y = x-Xh/x— \x\/x, 
where x^ is the scaled radius of the sphere of influence of a 
SMBH (see Appendix A for details). 

When the quasars at the mean redshift < z s >= 1-27 are 
lensed by foreground CDM halos of galaxies and clusters of 
galaxies, the lensing probability with image separations larger 
than A6* and flux density ratio less than q r is (Schneider et al. 
IT992T) 

/>(> A9, < q r ) = f dDh<KZ) dz f n(M,z)a-(M,z)B(M,z)dM, (7) 
Jo dz Jo 

where Z)l(z) is the proper distance from the observer to the 
lens located at redshift z. h(M, z) is the physical number density 
of virialized dark halos of masses between M and M + dM 
at redshift z given by Jenkins et al. (2001). The cross section 
cr(M, z) is mass and redshift dependent, and is sensitive to the 
flux density ratio of multiple images q T for SIS halos, as shown 
below in detail. For the magnification bias B(M, z), we use the 
results given by Li & Ostriker d2002> for both SIS+SMBH and 
NFW lense. For the former, if the contribution from a SMBH is 
ignored, the magnification bias for SIS only will be recovered 
(which is fisrs ~ 4.76). 

The cross section for the cluster-size NFW lenses has been 
well studied (Li & Ostriker I2002i . The lensing equation is 
y = x - fisg(x) / x and the multiple images can be produced only 
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brighter images 
fainter images 



Fig. 1. Image magnifications. Solid line: the brighter images 
(fj.+(y), x > 0), dashed line: the fainter images (fi-(y), x < 0). 
The left vertical line (dot-dashed) indicates the source position 
y CI at which q r = 10, and the right vertical line (dotted) indi- 
cates the position for q r = 10 4 . 

if \y\ < y a , where y a - is the maximum value of y when x < 0, 
which is determined by dy/dx = and the cross section in the 
lens plane is simply cr(M,z) = Try^jf. We point out that the 
flux density ratio would also result in a reduction of the value 
of the cross section for NFW lenses if we apply it in our cal- 
culations. However, since we model the cluster-size lens halos 
as an NFW profile, which contributes to lensing probabilities 
mainly at larger image separations, where no confirmed lens- 
ing events are found (i.e. the null results of the JVAS/CLAS for 
6 < AO < 15 ), if we adopt the mentioned above usual form of 
an NFW lens cross section, the predicted lensing probabilities 
will not be severely affected at smaller image separations. 

For galaxy-size SIS+SMBH modeled lenses, two images 
will always be produced for all values of the source positions 
\y\, and larger |j>| will produce a higher flux density ratio q t . So a 
finite value of q x will definitely limit the source within a certain 
corresponding position (which is also denoted by y CI ) and hence 
reduce the value of the cross section. The flux density ratio q t 
for the two images is the ratio of the corresponding absolute 
values of magnifications (Schneider et al. 1992 Wu 1996 1, q r = 
where 



My) 



iy_dy\ 
\ x dx) x>0 



and //_(;y) = 



ydy\ 
x dx) x<0 



(8) 



Soy cr is determined numerically by \fi+(y CI )\ = q T \[i-(y C r)\. 

We plot the magnifications (each as a function of source 
position y) both for the brighter images and the fainter images 



in Fig.n with the source positions at which q x = 10 and q T = 
10 4 explicitly indicated, and we have used a typical value xj, = 
10~ 3 . The sensitivity of the source position to the flux density 
ratio is obvious. The cross section for images with a separation 
greater than AO and a flux density ratio less than q r in the lens 
plane is (Schneider et al. 1992 Chen l2003b> 



cr(M,z) 




in) 

for AO < A0 ; 

for A6» < AO < A0 ya 

for AO > A0 Vct , 



(9) 



where §(x) is a step function, and M m ; n is the minimum mass of 
halos above which lenses can produce images with separations 
greater than AO. From Ea. (IA.5> . an image separation for any y 
can be expressed as A0(y) = ^Ax(y) / Z)j\ where Ax(y) is the 
image separation in lens plane for a given y. So in Eq.(|9j, the 
source position y^g, at which a lens produces the image sep- 
aration AO, is the reverse of this expression. A0q = A0(0) is 
the separation of the two images which are just on the Einstein 
ring; A0 ya = A0(y cr ) is the upper-limit of the separation above 
which the flux ratio of the two images will be greater than 
q T . Note that since Mum (Mis) is related to AO through ifo in 
equation (|6) and cr^, = GMuM/2r v [ T , we can formally write 
Mum = MuM(A0(y)) and determine M m ; n for galaxy-size lenses 
by M min = M DM (A% cr )). 



4. Discussion and conclusions 

The lensing probabilities predicted by Eq. and calculated 
from the combined JVAS/CLASS data are compared in Fig. 13 
Four cases with different flux density ratios q T or cooling scale 
M c or with and without central SMBH are investigated to show 
their effects on lensing probabilities. In each case (see the cor- 
responding panel of Fig. 0, four different values of the cos- 
mic equation of state ai: a> = -1 (cosmological constant), 
<jj = -2/3, a> = —1/2 and u> = -1/3, are chosen to show the 
constraints from lensing probabilities. As in our previous work 
(Chen 2003b I, the observed lensing probabilities are calculated 
by P bs(> AO) = N(> A<9)/8958, where N(> AO) is the number 
of lenses with separation greater than AO in 13 lenses. 

When flux density q r = 10 (as allowed by the JVAS/CLASS 
well-defined sample) and M c = 10 13 /i _1 Mq, one can see from 
the panel (a) of Fig. |3 that none of the four values of to is 
able to predict enough lensing probabilities to match the ob- 
servations. Note that the value of the cooling mass scale used 
in this case is preferred by Kochanek & White (2001 ) and Li 
& Ostriker (2002i, and is also close to the cutoff mass of ha- 
los below which cooling of the corresponding baryonic com- 
ponent will lead to concentration of the baryons to the inner 
parts of the mass profile (Rees & Ostriker [19771 Blumenthal 
et al. 119861 Porciani &Madau l2000l i. So a larger value of M c 
(> 1O 13 /i~'M0) is needed for predictions to match observa- 
tions. When we attribute too few predicted probabilities in 
panel (a) of Fig. [2] to the small value of q r (= 10), the quite 
large value of q T (= 10 4 ) also fails to match the observations, 
although the increment of the probabilities is also obvious, as 
shown in panel (c). The reason is that we have used different 
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Fig. 2. The integral lensing probabilities with image separa- 
tions larger than AO and flux density ratio less than q r , for 
quasars at mean redshift < z s >= 1.27 lensed by NFW (Mdm > 
M ) and SIS+SMBH (M DM < M c ) halos. The typically se- 
lected values of the cosmic equation of state to for each curve, 
the flux density ratio q r and cooling scale M c are explicitly indi- 
cated. The observed lensing probabilities of JVAS/CLASS are 
plotted as a histogram in each panel. 



SIS lenses considerably without a small central point mass. On 
the other hand, our calculations shows that the central SMBH 
can be detected only if q, > 10 3 , this value is possible for some 
confirmed radio-loud lens systems (Rusin & Ma 2001 ), but no 
sample suitable for analysis of the lens statistics is available 
with such a high flux density ratio. The sensitivity of lensing 
probability to the flux density ratio is obvious when we com- 
pare panel (b) with panel (d) of Fig. [2] In the latter, q r = 10 4 , 
the effect of which is close to that when flux density ratio is as- 
sumed to be infinite (i.e., no constraints on the flux density ratio 
are taken into account, as in most previous work; see Fig. 0, 
and other parameters, including M c , are the same as panel (b). 
Obviously, for q T = 10 4 , the "right" cosmic equation of state 
to match observations is to — -1/2. The present time deceler- 
ation parameter for a flat, dark energy dominated Universe is 
qo = 1/2 + (3/2)<jOq, and the accelerating Universe requires 
qo < 0, which immediately gives to < -1/2.1 (we are using 
Q.q = 0.7). So the fit value of to = -1/2 in panel (d) is still 
within the range required by the accelerating Universe. Clearly, 
when flux density ratio q T changes from 10 to 10 4 , the corre- 
sponding fit value of the cosmic equation of state will change 
from co = — 1 to o = -1/2. 

In summary, we revisit the two populations of lens halo 
model with mass distribution NFW (Mdm > M c ) and 
SIS+SMBH (Mdm < M c ), to calculate lensing probabilities 
in flat, low-density cosmological models with different cosmic 
equations of state a>. The finite flux density ratio effect is signif- 
icant. A low value of q r can be compensated for by raising the 
cooling mass scale M c in fitting the predicted lensing probabil- 
ities to JVAS/CLASS observations. To determine the cosmic 
equation of state a>, the uncertainty in M c must be resolved. 
The contributions from galactic central SMBHs can be safely 
ignored. 



parameters and functions, e.g., the sensitive concentration pa- 
rameter c\. Although it is not explicitly displayed in the fig- 
ure, our calculations indeed show that if we set q x — 10 4 and 
M c = 1.5 x lO 13 /r'M0, the predicted lensing probabilities for 
co = -1 is able to match the observations quite well, which is in 
agreement with the result obtained by Sarbu et al. (2001 1, since 
we have used almost the same parameters and mass (and halo) 
functions. Note that when q r is quite large (> 10 3 ), the contri- 
bution from the central galactic SMBH cannot be ignored. In 
fact, it is just this contribution that compensates the finite value 
effect of q r - 10 4 , which make our results obtained in this case 
be the same as that of Sarbu et al. J2001> . whose flux density 
value is infinite. 

If we persist in the flux density ratio allowed by the 
JVAS/CLASS sample in our calculations, the cooling mass 
scale M c must have the value of 3 x 10 13 hr x Mq for the model 
with u — -1 to be able to match the observations, as shown in 
panel (b). Such a value of M c is close to that used by Porciani 
& Madau ( 2000 1, however, they used different parameters and, 
especially, they did not consider the finite flux density ratio ef- 
fect. Note that when q T < 10, the contributions from a galactic 
central SMBH can be ignored, no matter what the value of M c 
is, as shown in the top two panels of Fig.|2] This means that the 
finite small flux density ratio will reduce the cross section of 
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Appendix A: Lensing equation for SIS+SMBH 

In this Appendix, we give a detailed deduction of the lensing 
equation for SIS+SMBH model. Choose the length scales in 
the lens plane and the source plane to be 

^-HiJ^f and * = *>5p (A - 1} 

then the position vector in the lens plane and source plane is, 
respectively, £, = jc£o and r\ = yrj^. Note that we have used 
the different length scales between SIS and NFW lenses. The 
dimensionless surface mass density then is 

M.5\x) al 1 

Ec^g 2GX cr £ |*| 

So the lensing equation is 

y = x - m(x)/x, (A.3) 



6 



D.-M. Chen: Flux-limited strong gravitational lensing and dark energy 



with 



m(x) 



= 2 I K(x')x'dx 
Jo 



M. 



+ X. 



(A.4) 



y = \y\ and x = \x\, and their values will be extended to the 
negative for convenience in later calculations. We have used 
the known formula 6 2 (x) = 6(x)/2jtx to derive the results in Eq. 
iAA\ . In terms of the radius of the sphere of influence or Bondi 
accretion radius of a SMBH (e.g., Ferrarese & Merritt 2000 
Melia & Falcke l200ll Zhao et al. 12002b r\ = GM./cr 2 „ the 
lensing equation can be rewritten as 



Xh \x\ r h 
y = x with Xh = . 

X X 7T^0 



(A.5) 



The mass of each galactic central SMBH M. is related to 
the mass of its host galaxy Mdm with the formula given by 
Ferrarese ( 2002 1 



M. 







2.78 x 10~ 4 Mj 5 57 , 



(A.6) 



here we have defined a reduced mass of the halos as M15 = 
^DM/(10 15 h^ 1 M o ). On the other hand, the mass of a SIS halo 
is defined the same way as for NFW: M DM = An J. v " p(r)r 2 dr. 
Since the mass density of a SIS halo is p(r) = cr 2 l2nGr 2 , 
the so-defined mass is related to the velocity dispersion by 
cr 2 v = GMuM/2r v i r . Generally, the virial radius r v ; r is cosmo- 
logical model-dependent according to its definition Mdm = 
4;7R) v j r /3(z)r) ir /3, where 6 v \ r is the well-known density contrast 
of the virialized dark matter halos with the familiar value 
<5 vir = I871- 2 w 178 for a standard cold dark matter (SCDM) cos- 
mology (Q m =1), and p(z) is the mean mass density of the uni- 
verse at redshift z. We use <5 v j r = 200 in our actual calculations, 
since Jenkins et al. (2001 1 specifically stated that their formula 
gives better fits to mass function f(M,z) with <5 v h- = 178 re- 
gardless of the cosmological model (see Sarbu et al. 2001). So 
the scaled radius of the sphere of influence of the SMBH Xh can 
be written explicitly 



x h = 5.18 x 10" 



-2/3 



(A.7) 



x[Q m (l+z) 3 +n Q (l+z) 3 < 1+ ">]~ 
where = D^D^ S /D^ and we have used (for a flat Universe) 
H\z) 



Pfe) = p(0)- 



= p(0)[n ra (l+z) 3 + f2 Q (l+z) 3(1+w) ]. 



(A.8) 
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